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There have been numerous studies addressing the turbulent diffusion in the 
stratosphere and mesosphere during the last two decades. The motivation for 
such studies was the need for an understanding of the thermal and constituent 
structure of the middle atmosphere. Observational estimates of the horizontal 
and/or vertical diffusion were obtained using chemical release, rocket vapor 
trail, aircraft, balloon, and radar techniques (ZIMMERMAN and CHAMPION, 1963; 
KOCHANSKI, 1964; REES et al., 1972; LILLY et al., 1974; ROSENBERG and DEWAN, 
1975; CADET, 1977; VINCENT and STUBBS, 1977; MANSON et al., 1979; SATO and 
WOODMAN, 1982; FRITTS, 1984). Typical_jstimates of vertical diffusion near the 
mesopause ranged from ^200 to 600 si s . During the same period, a number 
of theoretical studies were performed to infer the level of vertical diffusion 
needed to account for observed constituent profiles (JOHNSON and WILKINS, 1965; 
COLEGROVE et al., 1966; CUNNOLD et al., 1980; ALLEN et al., 19|1 ;_|PRUZESE et 
al., 1984). These estimates were for the most part near 100 in s . Thus, 
there appears to be a discrepancy between the level of vertical diffusion 
required for the dissipation of gravity wave and tidal motions on the one hand 
and for the maintenance of observed temperature and constituent profiles on the 
other. The purposes of the present contribution are to outline a possible 
explanation of this discrepancy and to suggest the measurements that may help, 
verify this explanation. 

The theory relies upon laboratory and atmospheric observations that suggest 
that the saturation of internal gravity waves via either convective or dynamical 
instabilities results in the generation of thin layers of turbulence which act 
to mix the local environment (CADET, 1977; KOOP, 1981; SATO and WOODMAN, 1982; 
WAND et al., 1983; PHILBRICK, private communication, 1983). Provided that these 
turbulent regions remain associated with the unstable portion of the wave field 
(as observed), it is then appropriate to examine the consequences of a turbulent 
diffusion that varies throughout the wave field in a systematic manner. Such 
an approach was taken in the study by FRITTS and DUNKERT0N (1984). The results 
of that study and their implications for atmospheric observations will be 
described here. 

The evolution of the wave and mean fields in the middle atmosphere are 
described by the nonlinear, viscous momentum and thermodynamic energy equations 
and the continuity equation. Using the continuity equation, the mean and 
perturbation potential temperature equations may be written 

®t + ( ' r "~e r) z " ^ z + V ; z (i) 

and 
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where subscripts denote differential and is the horizontal component of 
turbulent diffusion assumed proportional to v. Then, consistent with the 
laboratory observations by KOOP (1981), we assume a symmetric distribution of 
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turbulent diffusion with a maximum at the point at which the total potential 
temperature is minimum (negative for convective instability). The total 
potential tanperature my be written as 

8 =8 (1 + ctcos <j>) (3) 

z z 

where <j> is the phase function 

<}> = kx + mz - kct M 

and a is a measure of wave amplitude (u > 1 for convective instability). 

We now multiply (2) by 8 1 , average horizontally, and note that all triple 
correlations a re sm all by virtue of the assumed negative correlation of v 1 and 
9 Solving w' 8', assuming steady-state saturation, and substituting into 
(1) yields, after some manipulation, 

T = ~~ (v8 +2 THT ' + 0 ,2 + e x ' 2 ) . (5) 
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Alternatively, if we write 

v** f(cos ) (6) 

and define a positive coefficient 

8 = _ (7) 

then (5) becomes 
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This expression relates all of the wave and turbulence contributions to the mean 
turbulence flux of potential temperature. In particular, we note that both wave 
and turbulence fluxes include a countergradient diffusion term due to the 
occurrence of a maximum turbulent diffusion in regions of minimum 0^. These 
countergradient terms reduce the total down-gradient flux of heat (and 
constituents ). expected on the basis of a uniform turbulent diffusion. 


The magnitude of this effect can be illustrated by assuming a distribution 
for turbulent diffusion of the form 


V o ^ 
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for n > 0 so that 
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The magnitudes of the turbulence and total (wave plus turbulence) fluxes are 
shown for various values of n (and negligible v^) in Figure 1 with dashed and 
solid lines, respectively. Note that even for relatively small n (broad 
distributions of turbulent diffusion), the countergradient terms results in an 
appreciable reduction of the total flux expected for a uniform turbulent 
diffusion (with n «* 0). In the case where v and contribute equally to the 
down-gradient wave flux, it is possible to show that, while the reduction of the 
total flux may be significant, the total flux remains down-gradient for all 
a and n > 0 (FRITTS and DUNKERTON, 1984). Additionally, this theory predicts 
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Figure 1. Variation of the turbulent ( -) and net ( ) 

rate of change of mean potential temperature normalized 
by the turbulent transport of the mean for the case of a 
vertical turbulent diffusivity alone. Note that even 
small values of n cause significant reductions in the 
net flux. 


Frandtl numbers in the mean and perturbation equations much larger than unit^ 
for nonuniform turbulent diffusion, as anticipated by CHAO and ’SCHOEBERL 
(1984). 

These results imply that, while a large turbulent diffusion may be required 
for the dissipation of gravity wave and tidal motions, a substantially smaller 
effective diffusion may act to mix mean gradients of potential temperature and 
constituents. The reason is that intense turbulence (and diffusion) is expected 
in precisely those regions of the gravity-wave field in which the total 
gradients of potential temperature and constituents are small or opposite to 
those of the mean state. Thus, it may not be sufficient to infer a mean 
vertical diffusion by assuming an average level of turbulence applied to the 
mean state. It may be necessary, instead, to infer the effective vertical 
diffusion from observations of the intensity and distribution of turbulent 
layers as well as the thermal structure. Such measurements could resolve the 
discrepancy between the level of turbulence required to dissipate observed wave 
motions and that needed to explain middle atmosphere constituent distributions. 
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